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Abstract 

We calculate the quasinormal modes and associated frequencies of 
the Bahados, Zanelli and Teitelboim (BTZ) non-rotating black hole. 
This black hole lives in 2+1-dimensions in an asymptotically anti-de 
Sitter spacetime. We obtain exact results for the wavefunction and 
quasi normal frequencies of scalar, electromagnetic and Weyl (neu- 
trino) perturbations. 
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1 Introduction 



When one is describing the evolution of some conservative system, one often 
considers a small perturbation or a small departure from a known solution 
of the system, and one generally arrives at a wave equation describing it. 
For a system with no explicit time dependence, one finds the normal mode 
solutions of the wave equation, satisfying certain boundary conditions, and 
one can then specify completely the perturbation as a linear superposition of 
these normal modes. In this case the operator associated to the perturbation 
is self-adjoint, the frequencies are real and the modes are complete. 

However, when one deals with open dissipative systems, as it is the case 
in this paper, such an expansion is not possible. Instead of normal modes, 
one considers quasi normal modes (QNM) for which the frequencies are no 
longer pure real, signalling that the system is loosing energy. Although QNMs 
are in general not complete and therefore insufficient to fully describe the 
dynamics (see [0], |2| and references therein), they nevertheless dominate the 
signal during the intermediate stages of the perturbation, being therefore 
extremely important. 

QNMs of black holes were first numerically computed by Chandrasekhar 
and Detweiler 0, and subsequent numerical simulations 0, |], |J showed that 
the amplitude is dominated, at intermediate times, by a ringing signal due 
to the QNMs. Aside from the pure mathematical interest, black hole's QNM 
calculations have been a very active field, and new methods, both numerical 
and analytical have been developed (see for a review). 

Up until very recently, all these works dealt with asymptotically flat 
spacetimes. In the past few years there has been a growing interest in asymp- 
totically AdS (anti-de Sitter) spacetimes. Indeed, the Banados-Teitelboim- 
Zanelli (BTZ) black hole in 2+1- dimensions M, as well as black holes in 
3+1 dimensional AdS spacetimes with nontrivial topology (see, e.g. |10|), 
share with asymptotically flat spacetimes the common property of both hav- 
ing well defined charges at infinity, such as mass, angular momentum and 
electromagnetic charges, which makes them a good testing ground when one 
wants to go beyond asymptotic flatness. Another very interesting aspect 
of these black hole solutions is related to the AdS/CFT (Conformal Field 
Theory) conjecture 0. For instance, due to this AdS/CFT duality, quasi- 
normal frequencies in the BTZ black hole spacetime yield a prediction for the 
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thermalization timescale in the dual two-dimensional CFT, which otherwise 
would be very difficult to compute directly. If one has, e.g., a 10-dimensional 
type IIB supergravity, compactified into a BTZ x S 3 x T 4 spacetime, the 
scalar field used to perturb the BTZ black hole, can be seen as a type IIB 
dilaton which couples to a CFT field operator O. Now, the BTZ in the 
bulk corresponds to to a thermal sate in the boundary CFT, and thus the 
bulk scalar perturbation corresponds to a thermal perturbation with nonzero 
< O > in the CFT. 

There has been some recent work on perturbations of Schwarszchild AdS 
spacetimes: Horowitz and Hubeny [|TTJ studied the QNM frequencies for 



scalar perturbations in 4, 5 and 7 dimensions. Wang et al [0, [13| studied 
scalar perturbations and QNMs on a Reissner-Nordstrom geometry, Chan 
and Mann |14]] studied the QNM frequencies for a conformally coupled scalar 
field. For work on BTZ black holes such as entropy of scalar fields, see []T5| 
and references therein. 

In this paper we shall consider the QNMs of the 3D non-rotating BTZ 
black hole ||. The non-rotating BTZ black hole metric for a spacetime with 
negative cosmological constant, A = — is given by 

2 2 

ds 2 = (-M + T -)dt 2 - (-M + ^dr 2 - r 2 d<j) 2 , (1) 

where M is the black hole mass. The horizon radius is given by r+ = M l / 2 l. 
We shall in what follows suppose that the scalar, electromagnetic and Weyl 
(neutrino) fields are a perturbation, i.e., they propagate in a spacetime with 
a BTZ metric. We will find that all these fields obey a wave equation and 
the associated QNM are exactly soluble yielding certain hypergeometric func- 
tions. As for the frequencies one has exact and explicit results for scalar and 
electromagnetic perturbations and numerical results for Weyl perturbations. 
To our knowledge, this is the first exact solution of QNMs for a specific model 
(see!). 

In section 2 we give the wave equation for scalar and electromagnetic per- 
turbations, and find the QNMs themselves and their frequencies. In section 
3 we find the wave equation for Dirac and Weyl (neutrino) perturbations and 
analyze their QNMs. 
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2 Perturbing a black hole with scalar and 
electromagnetic fields 

2.1 The wave equation 

In this subsection we shall analyze the scalar and electromagnetic perturba- 
tions, which as we shall see yield the same effective potential, and thus the 
same wave equation. 

First, for scalar perturbations, we are interested in solutions to the min- 
imally coupled scalar wave equation 

Vw = 0, (2) 

where, a comma stands for ordinary derivative and a semi-colon stands for 
covariant derivative. We make the following ansatz for the field $ 

^ = ^f(r)e~ iuJt e im ^ (3) 



where m is the angular quantum number. It is useful to use the tortoise 



coordinate r* defined by the equation dr* = — ^—pr, and given implicitly by 



r = -M 1 ' 2 coth(M 1/2 r,) , (4) 

with r* e ] — oo, 0], (r* = — oo corresponds to r = r + , and r* = corresponds 
to r = oo). 

With the ansatz (|3|) and the tortoise coordinate r*, equation @ is given 

by, 

d 2 f(r) 
dr* 2 

where, 



+ (u-V(r))f(r) = 0, (5) 



3r 2 M M 2 m 2 Mm 2 , . 

and it is implicit that r = r(r#). The rescaling to the the radial coordinate 
f = j and to the frequency Cj = ul is equivalent to take I = 1 in (||) and @, 
i.e., through this rescaling one measures the frequency and other quantities 
in terms of the AdS lengthscale I. 
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Now, the electromagnetic perturbations are governed by Maxwell's equa- 
tions 



F" v . v = ,withF^ = A 



A 



(7) 

where F^ v is the Maxwell tensor and is the electromagnetic potential. As 
the background is circularly symmetric, it would be advisable to expand A^ 
in 3-dimensional vector spherical harmonics (see [1£§ and [|T7|): 



g m (t,r) 

h m (t,r) 
k m (t, r) 



z im<t> 



(8) 



where m is again our angular quantum number, and this decomposition is 
similar to the one in eigenfunctions of the total angular momentum in flat 
space 
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However, going through the same steps one finds that the equation for 
electromagnetic perturbations is the same as the one for scalars, equation 
(§|). The reason is that in three dimensions the 2-form Maxwell field F = 



F^dx^ A dx v is dual to a 1-form d<&. 



2.2 QNMs for scalar and electromagnetic perturba- 
tions 

Although a precise mathematical definition for a QNM can be given, as a 
pole in the Green's function j7|, we shall follow a more phenomenological 
point of view. A QNM describes the decay of the field in question. For 
the equation (|5]) it is defined as a corresponding solution which (i) near the 
horizon is purely ingoing, ~ e tujr * , corresponding to the existence of a black 
hole, and (ii) near infinity is purely outgoing, ~ e~ lujr *, (no initial incoming 
wave from infinity is allowed). One can see that the potential V(r) diverges at 
infinity, so we require that the perturbation vanishes there (note that r = oo 
corresponds to a finite value of r*, namely r* = 0). This vanishing of the 
solution at oo will only be possible for a discrete set of complex frequencies 
to called quasinormal frequencies. 
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2.2.1 Exact calculation 



Puting I = 1 and using the coordinate r*, the wave equation 
form 



takes the 



lo — 



3M 



4sinh(M 1 /2 r *) 2 + 4cosh(Af 1 /2 r ») 2 " r cosh(M 1 /2 r< ,)2 



d 2 a{r) 
M 



+ 



a(r) = 



(9) 



On going to a new variable x 
be written as 



coehCMVar.) 2 ' X 6 t 0, ^ eC l Uation (D Can als ° 

4x(l - x)d 2 x a + (4 - 6z)0 x ^ + V(x)a = , (10) 



where 
V(x) 



1 



4x(l — rr) 



4a; 2 (1 



Af 



3x — z(l 



4m 2 x(l — x) 
M 



By changing to a new wavefunction y (see [JT8 ] for details), through 

(x - I) 3 / 4 



■0 



X 2M 



equation (|10| ) can be put in the canonical form []18 , |l£ 



x(l - x)y" +[c- (a + b+ l)x\y' - aby = 



with a = 1 + 



1 



and c 



(12) 



(13) 



which is a 



217172 2*7172'" x 2MV2 2M1/2 ' aliu L x A7i72' 

standard hypergeometric equation. The hypergeometric equation has three 
regular singular points at x = 0,x = l,x = 00, and has two independent 
solutions in the neighbourhood of each singular point. We are interested 
in solutions of (0) in the range [0,1], satisfying the boundary conditions of 
ingoing waves near x — 0, and zero at x — 1. One solution may be taken to 
be 

(14) 



y=(l- x) c ' a ' b F(c -a,c-b,c;x), 



where F is a standard hypergeometric function of the second kind. Imposing 



y — at x — 1, and recalling that F(a, b, c, 1) = S^^J 

a = —n , or b = —n , 



, we get 



(15) 
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with n = 0, 1, 2, ... , so that the quasi normal frequencies are given by 

uj = ±m-2iM 1/2 (n + 1). (16) 

The lowest frequencies, namely those with n = and m = had already 
been obtained by j2(J and agree with our results. 

2.2.2 Numerical calculation of the frequencies 

In order to check our results, we have also computed numerically the fre- 
quencies. By going to a new variable z = ~, h = ^- one can put the wave 



equation fl5|) in the form (see JTTJ for further details) 



S ^ Z) Jz^ Q + ^Tz® + U ^ )Q = ° ' (17) 

where 9 = e iur *a{r), s(z) = z 2 - Mz 4 , t(z) = 2Mz 3 - 2iuz 2 and u(z) = 
zrjtfl i ) w hh V given by @. Now, z e [0,h] and one sees that in this 

range, the differential equation has only regular singularities at z = and 
z — h, so it has by Fuchs theorem a polynomial solution. We can now 
use Frobenius method (see for example ||21|| ) and look for a solution of the 



form Q(z) = J2^Lo^n(uj)( z — h) n (z — h) a , where a is to be taken from the 
boundary conditions. Using the boundary condition of only ingoing waves 



at the horizon, one sees that a = 0. So the final outcome is that G can 
be expanded as 

oo 

e(z) = J £e n{u) (z-h) n . (is) 

n=0 

Imposing now the second boundary condition, = at infinity (z = 0) one 

gets 

oo 

£0 n(w) HO n = o. (is) 

71=0 

The problem is reduced to that of finding a numerical solution of the poly- 
nomial equation (|T9|). The numerical roots for w of equation ([19]) can be 
evaluated resorting to numerical computation. Obviously, one cannot de- 
termine the full sum in expression (]T9"|), so we have to determine a partial 



7 



sum from to N, say, and find the roots uj of the resulting polynomial ex- 
pression. We then move onto the next term N+l and determine the roots. 
If the method is reliable, the roots should converge. We have stopped our 
search when a 3 decimal digit precision was achieved. We have computed 
the lowest frequencies for some parameters of the angular quantum number 
m and horizon radius r+. The frequency is written as uj = uj r + iuji, where 
uj r is the real part of the frequency and Ui is its imaginary part. 

In tables 1 and 2 we list the numerical values of the lowest QNM frequen- 
cies, for m = and m — 1, respectively, and for selected values of the black 
hole mass. 



m = 




Numerical 


Exact 




UJ r 


-UJi 


UJ r 


-UJi 


l 

2 


0.000 


1.000 





1 


1 


0.000 


2.000 





2 


5 


0.000 


10.000 





10 


10 


0.000 


20.000 





20 


50 


0.000 


100.000 





100 


100 


0.000 


200.000 





200 


1000 


0.000 


2000.000 





2000 



Table 1. Lowest (n = 0) QNM frequencies for m — 0. 



m—1 




Numerical 


Exact 




uj r 


-UJi 


UJ r 


-UJi 


l 

2 


1.000 


1.000 


1 


1 


1 


1.000 


2.000 


1 


2 


5 


1.000 


10.000 


1 


10 


10 


1.000 


20.000 


1 


20 


50 


1.000 


100.000 


1 


100 


100 


1.000 


200.000 


1 


200 


1000 


1.000 


2000.000 


1 


2000 



Table 2. Lowest (n = 0) QNM frequencies for m—1. 
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The numerical results agree perfectly with ([16]), and one sees that the 
imaginary part of the frequency scales with the horizon whereas the real 
part depends only on the angular index m. 



3 Perturbing a black hole with Dirac and 
Weyl spinor fields 

3.1 The wave equation 

We shall develop Dirac's equation for a massive spinor, and then specialize 
to the massless case. The two component massive spinor field with mass 
Us obeys the covariant Dirac equation 

i^V^ -fi a V = Q, (20) 

where is the spinor covariant derivative defined by = <9 M + |^C7[a7&] > 
and uj^ b is the spin connection, which may be given in terms of the tryad e£. 

As is well known there are two inequivalent two dimensional irreducible 
representations of the 7 matrices in three spacetime dimensions. The first 
may be taken to be 7 = ia 2 ,7* = a 1 , and 7 2 = a 3 , where the matrices a k 
are the Pauli matrices. The second representation is given in terms of the first 
by a minus sign in front of the Pauli matrices. From equation (|20|), one sees 
that a Dirac particle with mass /j, s in the first representation is equivalent to a 
Dirac particle with mass — fi s in the second representation. To be definitive, 
we will use the first representation, but the results can be interchanged to 
the second one, by substituting /j, s — > —[i s . For Weyl particles, fi s = 0, both 
representations yield the same results. 

Again, one can separate variables by setting 



*i(t,r) 



e -^V m0 . (21) 



On substituting this decomposition into Dirac's equation (^) we obtain 

" m 2A jf ) r ^ + * AV2 ^ 2 + = (m + , (22) 

i(M — 2r 2 ) r 2 uj 
~ V 2A1/2 + iA 1/2 d r ^i + = (m + a*.)* 2 , (23) 
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where we have put A = —Mr 2 + |r, we have restored the AdS lengthscale 
/, and in general we follow Chandrasekhar's notation [22|. Defining f?i, R 2 , 
and rh through the relations 

1 = iA' 1/4 Ri , (24) 
^ 2 = A _1/4 i? 2 , (25) 
m = im , (26) 



we obtain, 



iA 1 / 2 

(d r , - iu)R 2 = — 5— (m - in s r)Ri , (27) 
iA 1 ' 2 

(d r , + iu))R! = — ^-(m + in,r)Ri . (28) 



Defining now u, Tj T 2 , and f* through the relations 



we get 



1/ = arctan(^) , (29) 
m 

R 1 = e^Ti , (30) 

R 2 = e^T 2l (31) 

= t* + — arctan(— M , (32) 
2uj m 

(df>, - iu)T 2 = WT 1 , (33) 

{d K - ilo)T 2 = WT 2 , (34) 



where, 



iA 1/2 (m 2 + u 2 r 2 ) 3/2 

r 2 (m 2 + tfx 2 ) + ^ V ; 

Finally, putting Z± = Ti ± T 2 we have 

(a^ + w 2 )z ± = i4^±, (36) 

with 

v± = jy 2 ±^. (37) 
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We shall be concerned with massless spinors (fi s = 0) for which f* = r*, and 
W = Thus, 

V^C;-M)±^C;-MV». (38) 

In the form ([57]) one immediatly recognizes that the two potentials V+ and 
VI should yield the same spectrum. In fact they are, in SUSY language, 
superpartner potentials derived from a superpotential W (see Once 
again, we can rescale r and take / = 1, by measuring everything in terms of 
1. 



3.2 QNMs for Weyl perturbations 

Similarly, the wave equation (|36D for Weyl (until recently also called neutrino) 
perturbations may be put in the form 



d,Z ± + 



uj 2 — m 



m ±Ml/2 sinh(M 1 / 2 r,) 2 ' 



cosh(M 1 /V,) 2 cosh(M 1 /2 r ,) 2 



Z ± = . (39) 



Going to a new independent variable, x = — sinh(M 1//2 r^), x e [oo, 0], we can 
write 



J 2 (l+X 2 ) _ m 2 J_ ;n r 



in I 

(l + x 2 )Z" + xZ'+[ M 1+ ^ 2 M1/2 }Z = 0. (40) 
By changing the wavefunction Z to x 

( M!/2 2m _ rctan 

X = e 2M 1 / 2 21 w ; (41) 

we have 

(1 + x 2 ) X " + (^ + x) X ' + (^)x = . (42) 

On putting s = s e [|, zoo], we have again the hypergeometric equation 
(|T3"D, with a = ]^T72, & = — jjrj2, an d c = ~ ± ^^72, so that the solution to 
the wave equation is again specified around each singular point, and is given 
by the analytic continuation of the standard hypergeometric function to the 



complex plane [|T^, [IIJ . 



Since infinity is located at s — 5, there is no easy way to determine the 
QNM frequencies, so we have to resort to numerical calculations. If we put 
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([55]) in the form (O) one again sees that it has no essential singularities, so 
the numerical method just outlined in the previous section may be applied. 
Moreover, since V + and V_ have the same spectrum GBfl and the same QNM 



frequencies |22| we need only to workout the frequencies for one of them. In 
table 3 we present the numerical results for the QNM frequencies for neutrino 
perturbations and for selected values of the black hole mass. 



m = 1 




Numerical 


M 1 / 2 




—U>i 


2 


0.378 


2.174 


5 


0.316 


5.027 


10 


0.224 


10.006 


50 


0.099 


50.001 


100 


0.071 


100.000 


500 


0.0316 


500.000 



Table 3. Lowest QNM frequencies for m — 1. 

For large black holes one can see that the imaginary part of the frequencies 
scale with the horizon (r + = M 1 / 2 ), just as in the scalar and electromagnetic 
case. We have also computed some higher modes, and the real part of the 
frequency u r , does not seem to depend on which mode we are dealing with, 
just as in the scalar and electromagnetic case. 



4 Conclusions 

We have computed the scalar, electromagnetic and neutrino QNM of BTZ 
black holes. These modes dictate the late time behaviour of the fields. In all 
cases, these modes scale with the horizon radius, at least for large black holes 
and, since the decay of the perturbation has a timescale r = — , this means 
that the greater the mass, the less time it takes to approach equilibrium. 
We have also found that for large black holes, the QNM frequencies are 
proportional to the black hole radius. Since the temperature of a BTZ black 
is proportional to the black hole radius, the QNM frequencies scale with the 
temperature, as a simple argument indicates fllfl . For the study of QNM 
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of 3+1-dimensional spherical, 
Lemos 



as well as the toroidal black holes found by 
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sec 
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Is there, for small black holes, any relation with these quasinormal modes 
and critical phenomena as speculated by Horowitz and Hubeny |TTJ ? Though 
that would be an extremely interesting relation, we still are not able to answer 
that. 
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